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Abstract:  ^  initio  calculations  of  isomeric  carbon  clusters  C„° 
and  Cn^,  n  =  2-4,  yield  structures  and  energies  similar  to  previous 
reports,  although  five  (not  two)  C**  structures  have  local  energy 
ninima.  Dissociative  ionization  of  structurally  varied  precursors 
/as  used  to  prepare  C3  and  C^  ionic  and  neutral  isomers;  however, 
•cheir  mass  spectra  from  collisionally  activated  dissociation  (CAD) 
and  neutralization-reionization  (NR)  under  a  wide  variety  of 
conditions  are  indistinguishable,  indicating  only  one  isomer  or  the 
same  mixture  of  isomers.  Likewise,  CAD  and  NR  spectra  of  C«*  and  0^° 
from  ^®CH2=CHCH=^®CH2  and  Cg^  and  €3°  from  CHg-^^CHCHg  show  complete 
13c/ 12c  scrambling.  CAD  cross  sections  are  consistent  with  C/  - 
ions  as  mainly  linear  isomers  and  Cj*  ions  from  cyclic  precursors  as 
mainly  cyclic.  Product  abundances  from  the  unimolecular 
dissociation  of  Cn",  C„' ,  and  C/  allow  the  selection  of 
thermodynamic  data  that  should  be  of  higher  relative  accuracy,  such 
as  11.4  eV  for  the  C3  ionization  energy  from  reported  values  of 


10.0-13.0  eV. 
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Carbon  clusters  were  first  studied  spectroscopically  using 
emission  from  interstellar  clouds  and  other  cosmic  environments^-^ 
and  by  mass  spectrometry  of  carbon  vapors.^  As  recently 
reviewed/’^  the  chemistry  of  molecules  and  ions  in  plasmas 
(carbon  arcs,  laser  ablation)^'®  and  in  the  formation  of 
polynuclear  aromatics,  diamond  films, ^  and  soot^*®  has  been  the 
subject  of  extensive  theoretical  and  experimental  research.  Of 
the  latter,  most  studies  have  measured  ionic  species;  in  this 
study  neutralization-reionization  mass  spectrometry  (NRMS)®  is 
used  to  investigate  the  isomers  and  energies  of  C^®,  as  well  as 
Cn^  and  Cn',  clusters  for  n  <  7  by  mass  selection  and 
neutralization  of  Cn^  ions  from  linear,  branched,  and  cyclic 
molecules. 

Theoretical  calculations  hav~  predicted  for  small 
species  that  odd-numbered  clusters  have  lower  electron 
affinities,  closed  shell  ground  states,  and  higher  stabilities; 
for  the  cationic  counterparts,  a  similar  order  of  stability 
was  indicated,  whereas  for  Cn‘  anions  an  opposite  trend  was 
found. Although  stable  C^®  and  C/  (n  <  10)  structures  were 
first  predicted  to  be  non-cyclic,  more  recent  ab  initio 
calculations  show  stability  also  for  cyclic  isomers, with  the 
two  isomers  close  in  energy  for  and 

In  experimental  studies,  carbon  vaporization,  with  and 
without  postphotoionization,  yields  mass  spectra  with  more 
abundant  odd-numbered  Cn^  and  even-numbered  Cn*  clusters, '‘’^®*^®  as 
predicted  by  theory.  This  is  more  indicative  of  the  relative 
stabilities  of  Cn*  and  Cn*  than  of  Cn®;  as  pointed  out  for  Cn* 
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spectra  by  Bowers  and  co-workers,^®  the  actual  Cn°  abundances 
measured  by  mass  spectrometry  depend  significantly  on  ionization 
cross  sections  and  stabilities.  However,  the  high  relative 
stability  of  €3°  is  supported  by  its  preferential  loss  from  both 
4.16-18  Mass-selected  clusters  undergoing 

metastable,  collisionally  activated,  and  laser  dissociation'' 
also  favor  formation  of  Cn^  species  of  odd-numbered  n  values. 

Concerning  isomeric  characterization,  evidence  for  open 
chain  Cn°  and  has  come  from  photoelectron, infrared,^®  and 
microwave^*^^  spectroscopy,  low  energy  Cn''’  dissociations,^®  and  ion 
mobility  values. However,  stable  cyclic  structures  are 
indicated  for  neutrals  and  C^*  ions  by  coulomb  explosion 

techniques,^®  and,  from  the  same  laboratory,  for  €3“  -  €7°  neutrals 
by  electron  affinity  measurements these  values  are  much 
lower  than  those  reported  earlier,^®  possibly  because  the  cyclic 
isomers  are  present  in  <10%  concentration,”®  with  this  value 
dependent  on  laser/graphite  interaction  conditions.” 
Additionally,  cyclic  is  indicated  by  infrared  spectra  using 
an  argon  matrix,”  and  a  second,  presumably  cyclic,  isomer  by 
differences  in  bimolecular  reactivity”  and  ion  mobility.”  All 
of  these  studies  utilized  high -temperature  formation  from 
elemental  carbon  of  these  clusters,  whose  formation  entropy 
should  favor  linear  isomers.'* 

In  a  novel  approach  that  could  avoid  this  entropic 
preference,  Lifshitz  and  co-workers”  utilized  dissociative 
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electron  ionization  of  cyclic  perchlorinated  molecules  to 
generate  with  n  =  3, 5, 6,7  and  >10,  but  found  no  evidence  that 
this  method  gave  different  isomeric  structures  than  those 
generated  by  carbon  vaporization.  Here  we  have  extended  this 
approach  in  an  attempt  to  prepare  linear,  branched,  and  cyclic 
Cn*  ions;  €3^  and  have  been  studied  the  most  extensively 
because  their  isomers  are  predicted  to  be  stable  and  of  closely 
similar  energies. Further,  with  NRMS®-^®  these  mass- 
selected  cations  are  then  utilized  to  form  the  corresponding 
neutrals  and  anions.  Product  abundances  from  their  competitive 
unimolecular  dissociations  are  used  to  evaluate  basic 
thermodynamic  values  proposed  for  0^°,  and  (Table  I)  /'5.27,28 

for  which  major  disagreements  still  exist.  Additional 
theoretical  calculations  have  been  carried  out  to  guide  these 
experiments. 

Experimental  Section 

Using  a  tandem  double-focusing  (EB-EB;  E  =  electrostatic,  B 
=  magnetic  sector)  mass  spectrometer  described  in  detail 
elsewhere,^®  ions  formed  by  70-eV  electrons  are  accelerated  by  10 
keV,  mass  selected  by  MS-I  (EB) ,  and  undergo  charge  exchange 
collisions  with  a  neutral  target  gas  in  a  first  collision  cell 
(Cls-I)  to  form  fast  neutrals.  The  calcium  neutralization 
experiments  employed  a  special  furnace  for  vaporization 
temperatures  up  to  800®  C  and  special  care  to  minimize  Hg 
background  (a  reason  that  xenon  was  used  as  a  high  IE  target 
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instead  of  Hg)  Residual  ions  are  deflected  electrostatically, 
and  the  fast  neutrals  are  made  to  undergo  collisional 
reionization  in  Cls~II  to  yield  cations  or  anions;  negative  ions 
are  also  formed  by  charge  reversal  of  the  precursor  cations  with 
CgHg  in  Cls-II.  The  resulting  ions  are  mass  separated  in  E-II 
and  detected.  For  cross  section  values,  the  neutral  beam  flux  is 
measured  using  a  retractable  channeltron  multiplier  at  Cls-III 
before  E-II.  Following  proposed  conventions,^®  the  designation 

Na(85%T)/O2(70%T)  indicates  a  spectrum  from  cation 
neutralization  with  Na  at  85%  precursor  transmittance,  followed 
by  residual  ion  deflection  (slash)  and  reionization  to  cations 
with  O2  at  70%  transmittance.  Data  acquisition,  reduction,  and 
computer  control  employed  a  PC-based  computer  system.*®*^ 

Specific  ion  precursors  included:  hexabromobenzene  (£- 
CgBre) ,  hexachlorobenzene  (c-CgClg) ,  benzene  (c-CgHg) ,  hexachloro- 
cyclopentadiene  (c-CsClg)  ,  1,3-butadiene  (n-C^Hg) ,  hexachloro-1, 3- 
butadiene  (xi“C4Clg) ,  l,4-dibromo-2-butyne  (ri-C^H^Brz)  ,  l-bromo-2- 
methylpropene  (i-C4H7Br) ,  l,2-dibromo-2-methylpropane  (i-C*H8Br2) , 
propene  (n-CsHg) ,  and  cyclopropane  (c-CsHg) .  H2^®C=CHCH=^®CH2  was 
synthesized  adopting  the  procedure  for  the  corresponding  ^*C- 
isotopomer;®^  all  other  compounds,  including  CH2=^®CHCH3,  were 
obtained  commercially. 

Theory,  standard  ab  initio  molecular  orbital  calculations®^ 
were  carried  out  using  the  GAUSSIAN  88  program®®  with  geometry 
optimizations  at  the  HF/6-311G(d)  level  using  the  Quadratic 
Convergence  technique  (C2^  and  structures  were  also  optimized 
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at  the  UMP2/6-311 (d)  level).  Harmonic  vibrational  frequencies 
were  used  both  to  characterize  stationary  points  on  the  surface 
as  minima  (all  frequencies  real,  representing  equilibrium 
structures)  or  first-order  saddle  points  (one  imaginary 
frequency,  representing  transition  structures) ,  and  (after 
scaling  by  0.89)  to  calculate  zero-point  vibrational 
contributions  to  relative  energies.  Improved  relative  energies 
were  obtained  on  the  HF/6-311G(d)  optimized  geometries  using 
Moller-Plesset  perturbation  theory  (frozen  core)  terminated  at 
second  (MP2)  and  fourth  (MP4)  orders. 

Results  and  Discussion 

Theory.  Our  ab  initio  calculations  (Table  II)  were  designed 
to  search  for  new  stable  isomers  to  guide  the  experimental 
investigation.  A  broader  exploration  of  the  energy  surface  was 
sought  here  by  utilizing  a  larger  basis  set,  6-311G(d) .  Ab 
initio  calculations  for  carbon  clusters  are  notoriously  difficult 
and  often  require  multi-configuration  methods  in  order  to  obtain 
correct  electronic  states  and  isomer  relative  energies. 

However,  this  was  impractical  in  these  extensive  calculations,  so 
that  some  stable  isomers  and  electronic  states  were  not 
identified.  As  found  previously, triplet  excitations  in  the 
MP4  treatment  exaggerate  the  relative  stabilities  of  some  isomers 
(Table  II) ,  but  yield  heats  of  atomization  that  are  in  better 
agreement  with  the  experimental  data  and  with  previous 
calculations  that  include  triple  excitation. 

Repeated  optimizations  of  the  doublet  found  only  the 
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state,  while  the  multi-configuration  calculations^^  find  this  as 
an  excited  state  275  kJ  mol'^  above  the  first  doublet  state. 

Similarly,  our  C2*  quartet  state  (*ng) ,  bond  length  1.222  A, 
corresponds  to  the  first  excited  state,  116  kJ  mol'^  above  the 
previously  identified  X*2g*  ground  state. For  the  Cj®  neutral 
found  previously,  the  global  minimum  is  the  state  (bond 
length  1.24  A),  while  the  triplet  corresponds  to  the  c^S„*  excited 
state;  the  triplet-singlet  energy  difference  (101  kJ  mol'^)  at 
MP4(SDQ)  level  is  close  to  that  published  (113  kJ  mol‘^)  but 
the  lowest  triplet  state  (^IIu)  found  previously'^-^  was  not 
identified  in  our  UHF  geometry  optimization.  The  cyclic 
isomer  has  been  found  by  many  recent  multi-configuration 
calculations^*  to  be  8-25  kJ  mol'^  more  stable  than  the  linear 
We  only  find  here  the  latter,  a  bent  isomer,  and  a 
linear  quartet  state  as  local  minima.  Neutral  C3  calculations 
found  the  linear  singlet  bond  length  1.275  A,  and 

triplet  (^IIu)  as  ground  and  excited  states,  with  an  energy 
difference  of  204  kJ  mol'^,  MP4(SDTQ),  close  to  the  calculated 
(197  kJ  mol'^)^®  and  experimental  (203  kJ  mol'^)®^  values. 

Two  stable  structures  were  identified  previously  by 
theory;*'^^  here  an  additional  three  were  found  as  local  energy 
minima  (Table  II  and  Figure  1) .  Of  the  species  already 
identified, we  find  the  linear  doublet  to  be  the  most 

stable  C/-  species  at  the  MP4(SDQ)  level  of  theory,  while  the 
linear  quartet  3^  (*Bg)  of  a  slightly  distorted  zigzag  geometry 
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(DoOh  -  Cjh)  is  the  most  stable  isomer  at  MP4(SDTQ).  The  cyclic  '“A" 
quartet  2*  is  another  stable  structure  of  a  slightly  distorted 
rhombus  geometry  (Djh  -♦  C,)  .  These  distortions  are  likely 
artifacts  of  the  UHF  geometry  optimization.^'  The  calculated 
quartet  states  show  significant  spin  contamination  (Table  II) , 
suggesting  the  existence  of  other  closely  spaced  high 
multiplicity  states.  Two  other  Ct*-  species,  the  cyclic  trapezoid 
4*  state)  and  the  branched  isomer  5*  (^Bi  state)  ,  although 
significantly  less  stable,  also  represent  minima,  motivating  our 
experimental  search  for  further  stable  isomers. 

Our  calculations  of  the  C4  neutral  system  (Table  II)  found 
no  new  isomers.  Geometry  optimizations  led  to  the  linear 
triplet,  bond  lengths  1.297  (end)  and  1.274  A,  identified  as  the 
^Ilg  state,  but  not  to  the  corresponding  ground  state^^'^'  lying 
125-145  kJ  mol'^  below.  We  also  found  rhombic  (^Ag)  and  linear 
(Tg*)  singlet  isomers,  138  and  20  kJ  mol'^  more  stable  than  the 
triplet  (^Ilg)  ,  at  the  MP4(SDQ)  level  of  theory,  consistent  with 
previous  calculations.^'  However,  the  linear  singlet  energy  is 
too  low  in  the  MP4(SDTQ)  calculations  (Table  II). 

Isomer  Characterization  from  CAD  and  NR  Spectra. 

Collisionally  activated  dissociation,  CAD  [02(50%T)]  of  10  keV  C2* 
-  clusters  gave  the  Figure  2  spectra;  the  previously  reported 
CAD  (air)  spectra  of  8  keV  Cg'^,  ,  and  C7*  ions  are  in 

agreement. However,  our  and  C*'^  spectra  are  independent 
(±10%)  of  the  wide  variety  of  precursor  molecules  used.  Further, 
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these  spectra  were  not  significantly  changed  by  collisionally 
activating  [He(50%T) ,  first  field-free  region]  the  precursor  ions 
to  increase  their  average  internal  energy.  The  CAD  spectra  of 
Cg*  from  fi-CgClg  and  s-CgBrg,  and  of  c/  from  toluene  and 
cycloheptatriene  also  were  closely  similar. 

NR  mass  spectra  have  been  used  successfully  to  differentiate 
a  multiplicity  of  isomers  for  cations  whose  CAD  spectra  are 
closely  similar,  such  as  those  of  C^Hg  and  C^H*.^®  However,  NRMS 
studies  with  the  wide  variety  of  Cn^  precursors  (Figures  3,4) 
found  no  isomeric  differences  outside  experimental  error.  The 
charge  exchange  spectra  [*NR'  C6H6(30%T) ,  Figure  3]  gave 
similar  (±20%)  data  for  Cj^  and  C^*,  irrespective  of  the  precursor 
used,  and  this  was  also  true  (±10%)  for  neutralization  of  Cn^ 
with  Xe,  Na,  or  K  followed  by  reionization  [02(70%T),  Figure  4]. 

All  C/  spectra  from  xenon  neutralization  (70%T;  IE  =  12.1 
eV)  are  dominated  by  the  reionized  precursor  peak,  with  fragment 
peak  abundances  closely  similar  to  those  in  the  corresponding  CAD 
spectra;  this  is  consistent  with  the  formation  of  these  fragments 
by  dissociation  of  cations  after  reionization,  not  by  neutral 
dissociation.®®  More  excited  neutrals  can  be  formed  by  using  a 
target  of  lower  ionization  energy  (IE) ,  thus  removing  less  energy 
from  the  ion  in  the  neutralization  process.®-^®  Using  Na  or  K  (IE 
=  5.1  or  4.3  eV)  greatly  reduced  the  reionized  precursor  peak  for 
Ci  -  Cg  (Figure  4) ,  consistent  with  a  corresponding  extent  of 
neutral  dissociation,  but  the  fragment  abundances  were  still 
independent  (±10%)  of  the  precursor  molecule  used  for  all 
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ions  studied. 


The  surviving  precursor  concentration  from  an  intermediate 
extent  of  0*°  dissociation  should  yield  a  more  sensitive  test  of 
isomeric  identity  than  the  abundances  of  the  Cj.3°  products,  as 
dissociation  could  be  preceded  by  isomerization  of  the  C*° 
precursor.  Neutralization  with  calcium  (IE  =6.1  eV)  did  provide 
partial  dissociation  (Table  III) ,  but  spectra  from  five 

different  precursors  still  had  closely  similar  abundances  of 
reionized  C4*  as  well  as  C1.3*  fragment  ions. 

As  a  further  test  of  the  stability  of  the  linear  isomer  of 
C/  and  €4°,  [1,4--^^C2]-1, 3-butadiene  was  used  as  a  precursor  in  an 
attempt  to  generate  Its  CAD  [02(70%T)]  and  NR 

[K(90%T)/02(70%T)  ]  spectra  showed  intensity  ratios  for  m/S. 
24:25:26  of  1.3:4. 7:1  and  1.1:4. 3:1,  respectively,  and  ratios  for 
m/£  37:38  of  0.99:1  and  1.00:1,  respectively.  For  these  peak 
abundances,  single  bond  dissociation  of  without 

isomerization  would  yield  0:1:0  and  1:0,  respectively:  with 
complete  scrambling,  these  ratios  should  be  1:4:1  and  1:1. 
Although  isomerization  is  thus  nearly  complete,  this  does  not 
necessarily  indicate  that  the  C4*  ion  is  not  linear  or  that  the 
linear  ion  is  unstable,  as  isotopic  scrambling  could  occur 
during  formation  of  the  €4^  ion  by  dissociative  ionization  of  the 
butadiene.  Similarly,  the  CAD  [He(70%T)],  NR  [K(90%T)/02(70%T)  ] , 
NR  [K(90%T)/He(70%T) ] ,  and  charge  reversal  [C5H6(30%T)]  spectra 
of  propene-2-“C  also  showed  similar  scrambling.  The  intensity 
ratios  of  lo/z  24:25  in  these  spectra  were  found  to  be  1:2.1, 
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1:2.0,  1:2.2,  and  1:2.1,  respectively;  complete  scrambling  should 
give  1:2.  Again,  the  and/or  €3°  structures  are  not 
necessarily  non-linear  or  unstable.  However,  all  of  the  above 
data  are  consistent  with  the  formation  of  both  and  from 
the  various  precursors  as  single  isomers  or  similar  mixtures  of 
isomers. 

CAD  and  NR  Efficiencies.  As  reported  earlier , the 
relative  yields  of  in  CAD  spectra  (Figure  2)  are  consistent 
with  the  expected  higher  stability  of  the  odd-numbered  cationic 
clusters.  The  relative  efficiencies  for  CAD  dissociation  of  C^* 
(Table  IV)  also  reflect  this  stability  relationship,  assuming 
increasing  dissociation  cross  section  with  increasing  size,  as 
found  for  metastable  ion  (-10'®  s  lifetime)  dissociation. 

Although  mass  discrimination®®  could  affect  the  summed  product 
cation  abundances  from  CAD  of  Cj,*,  the  value  for  C/  is 
substantially  below  those  of  C^*,  and  Cg^.  Metastable  ion 

studies  found  a  small^^  or  significant^®  decrease  in  the  cross 
section  for  C7'^  versus  Cg'^,  but  still  greater  than  the  value  for 
C5*  (C^*  was  not  measured)  .  For  ion-molecule  reactions  of  laser¬ 
generated  Cy'*’  ions,  two-thirds  were  found  to  be  unreactive  with 
D2,  indicating  that  the  majority  of  ions  have  the  cyclic 
structure.®®  The  relatively  low  degree  of  fragmentation  for  C7^ 

C<7'^  (Table  IV)  indicates  that  Cj*  is  unusually  stable  and/or  of 
smaller  physical  cross  section  for  undergoing  collision;  either 
is  consistent  with  a  possibly  higher  fraction  of  cyclic  isomers 
for  C7‘^  prepared  by  dissociative  ionization  of  cyclic  precursors. 
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As  a  corollary,  this  also  implies  that  Cj^  and  Cg^  are  mainly  the 
linear  isomers.  For  the  even-numbered,  less  stable  and  €5+ 
clusters,  the  increase  in  -*  C</  cross  section  values  should  be 
due  primarily  to  increasing  size,  consistent  also  with  a  linear 
structure  for  Ci,*. 

The  Cj,*  -*  C^°  neutralization  efficiency  (Table  IV)  decreases 

quite  regularly  from  Ci  (15.4%)  to  Cg  (4.9%);  the  increasing 

physical  cross  section  apparently  is  more  effective  in  producing 

scattering  than  electron  transfer.  The  0^°  -♦  C<„^  reionization 

efficiency  decreases  even  more  dramatically  from  (5.2%)  to  Cg 

(0.3%)  despite  decreasing  IE  values;  neither  the  C„^  ->■  C„°  nor 

-*  C„*  efficiencies  reflect  the  odd/ even  nature  of  the  cluster, 

consistent  with  the  tact  that  odd  clusters  are  more  stable  for 

both  Cn°  and  Cn"^.  In  contrast,  the  efficiency  for  forming  C„'  from 
.  stub/f/ty 

C^^as  well  as  the  Cn'^shown  by  the  values,  favors  the 

even-numbered  clusters,  as  expected.  Although  the  efficiency  for 
-♦  C„*  drops  substantially  from  C2  (0.90%)  to  Cg  (0.14%),  most 
of  this  appears  to  be  due  to  the  drop  in  neutralization 
efficiency,  C^^  -♦  C„°;  the  negative  effect  of  size  is  far  less  for 
Cn°  -♦  Cn'  than  for  the  much  higher  energy  C„°  C^^  process. 

Cn“  Heats  of  Formation  (aH^).  The  aHj  values  of  C^®  from 
experimental  and  theoretical  studies  (Table  I)  were  checked  for 
their  ability  to  predict  competitive  product  abundances  from 
unimolecular  dissociations  (Figure  4)  after  neutralization.  This 
assumes  minimal  entropy  requirements;  energy  release  measurements 
for  Cn"*^  metastable  ion  dissociations  indicate  negligible  reverse 
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activation  energy.^®  ”  Comparisons  can  only  be  qualitative 
because  of  serious  mass  discrimination  against  smaller  product 
ions.  For  example  (Figure  4)  [Cj^] ,  =  717  kJ  mol'^,  in  the 

NR  spectrum  should  even  be  higher  than  [C2*”] ;  the  smaller  C„° 
products  have  fewer  vibrational  degrees  of  freedom  (none  for 
Ci°)  ,  so  that  a  higher  proportion  of  excess  excitation  energy 
could  become  translational  energy,  increasing  loss  of  the  smaller 
product.  Offsetting  this  somewhat,  these  do  have  higher 
reionization  efficiencies  (Table  IV) ,  and  their  abundances  are 
also  increased  by  secondary  product  dissociation,  as  shown  by  K 
versus  Na  neutralization. 

For  NR  data  (Figure  4  and  Table  III) ,  formation  of  + 

C3  is  somewhat  favored  over  that  of  2  C2.  Using  AHf(Ci°)  =  717  kJ 
mol'^  as  the  reference  value, the  pairs  of  AHf(C2°,  €3°)  values 
(Table  I)  of  (832,  837), (815,  801),““  and  (782,  773)  predict 
that  Cj®  formation  is  lio,  115,  and  74  kJ  mol'^,  respectively, 
more  endothermic  than  that  of  Ci°  +  0^°;  even  the  74  kJ  mol'^  could 
be  high.  For  Cj®  dissociation,  using  AHf(C2°,  Cg®)  of  782  and  773 
with  AHf(C4®)  =  971,  1022,  and  1012  kJ  mol"^  predicts 
endothermicities  of  133,  184,  and  174  kJ  mol"^  favoring  the 
formation  of  Cj®  +  C3®  over  that  of  Cg®  +  again  the  Figure  4 

data  would  appear  to  favor  the  lowest  value  (the  C^,®  reionization 
efficiency.  Table  IV,  is  half  that  of  Cj®)  .  For  Cg®  dissociation 
using  aH£(C2®,  C3®,  C^®)  of  782,  773,  and  971  kJ  mol"^  predicts  an 
endothermic ity  of  207  kJ  mol'^  favoring  formation  of  2  C3®  over 
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that  of  C2°  +  Ci,°;  this  value  appears  more  consistent  with  the 
Figure  4  data  than  those  of  258  and  248  kJ  mol'^  derived  with 
AHf(C^°)  =  1022  and  1012  kJ  mol'^.  Using  the  AHf(Cy)  values  of  979 
and  1030  kJ  mol‘^  with  these  aI^(C2°,  C^°)  values  predicts 
endothermicities  of  57  and  6  kJ  mol'^  favoring  formation  of  + 
C5°  over  that  of  +  C°‘,  to  agree  with  the  Figure  4  data  and 
subsequent  correlations,  979  kJ  mol‘^  has  been  chosen.  These 
values,  underlined  in  Table  I,  should  be  of  higher  relative 
accuracy  than  of  absolute  accuracy.  These  predict  minimum 
energies  of  dissociation  of  €2°,  652;  Cj®,  726;  0^°,  519;  C°,  576; 
and  Cg®,  343  kJ  mol"^,  consistent  with  higher  stability  for  the 
odd-n-value  neutral  clusters. 

0^.°  Ionization  Energies  (lE^  .  The  agreement  in  reported  IE 
values  (Table  I)  is  much  poorer  than  even  the  accuracy  presumed 
for  our  underlined  aHf  values.  For  example,  the  10.4  to  13.0  eV 
values  reported  for  IE  (€3°)  are  both  substantially  above,  and 
substantially  below,  those  of  IE(C2°)  and  IE(C4,°)  (Table  I).  For 
this,  comparison  of  the  relative  product  yields  in  the  NR  and  CAD 
spectra  should  be  useful,  as  these  should  be  affected  mainly  by 
the  relative  IE  values,  with  mass  discrimination  at  least 
qualitatively  similar.  Thus  for  the  CAD  spectrum  of  Cj^  -►  Ci°  + 

C2‘^  vs.  C2°  +  Ci^,  the  0.8  eV  predicted  lower  IE  of  C^”  should 
correspondingly  favor  C^*  formation  versus  that  of  C2*,  consistent 
with  the  CAD  (Figure  2)  versus  NR  (Figure  4)  spectra  of  C3^. 
Comparing  the  CAD  and  NR  spectra  of  C^"^,  [C3'^]/ [Cj"^]  and  [C3°]/[Ci°] 
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are  similar,  suggesting  comparable  IE(C3°,Ci°)  values.  The 
[C3*]/[C2*]  yield  is  several  times  that  of  [C3°]/[C2°],  again  in 
agreement  with  the  lower  IE (€3°)  values  found  by  theory,  but  not 
by  other  experiments  (Table  I) ;  this  is  also  borne  out  by  the 
C3/C2  values  in  the  CAD  and  NR  spectra  of  C5*.  The  C5*  spectra 
also  show  that  the  [C/]/[C3^]  yield  is  somewhat  greater  than  that 
of  tC4°]/[C3°],  indicating  IE(C4°)  <  IE(C3°)  ;  also  consistent  with 
this  are  the  data  of  the  CAD  and  NR  spectra.  For  the  Cg"^ 
spectra  the  [C5^]/[C/]  yield  is  somewhat  greater  than  that  of 
[^5°]/ [C<,°] ,  although  this  is  offset  by  the  poorer  reionization 
efficiency  of  C^°  vs.  C<,°  (Table  IV).  Thus  the  IE  values  of  C^® 
and  Cg®  should  be  comparable,  with  C2°  higher  and  and  lower, 
consistent  with  the  underlined  values  of  Table  I. 

Values.  The  heat  of  the  formation  values  of  the  Cp"^ 
species^  which  represent  the  sum  of  the  and  IE  values  of  the 
respective  neutrals,  can  be  used  to  predict  the  enthalpy 
differences  in  C^^  dissociation  products  of  the  CAD  spectra 
(Figure  2) .  The  best  agreement  with  these  data,  again 
recognizing  mass  discrimination,  was  obtained  using  the 
underlined  aH£(C„'^)  values  of  Table  I.  From  formation  of  C2^ 
is  predicted  to  be  less  favored  by  82  kJ  mol'^,  but  mass 
discrimination  should  greatly  reduce  [Cj^] ,  which  is  still  an 
abundant  peak  in  the  Cj^  CAD  spectrum.  From  Ci^  formation  is 
13  kJ  mol‘^  favored  over  that  of  C3*,  by  far  the  most  abundant 
peak,  with  C2^  less  favored  by  156  kJ  mol'^.  From  , 
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differential  enthalpy  values  are  C**,  67;  C3*,  0;  C2* ,  69,  and  Cj*, 
120  kJ  mol'^;  from  ,  these  values  are  C5*,  84;  C/,  141;  €3^,  0; 

C2  ,  276;  and  137  kJ  mol'^;  and  from  Cj* ,  these  values  are  Cg^, 

89;  C5*,  17;  C*'",  0;  ,  66;  Cj*,  152;  and  C*,  229  kJ  mol'^  These 

values  are  qualitatively  (±30  kJ  mol‘^)  consistent  with  the 
Figure  2  spectra.  These  predict  minimum  energies  of  dissociation 
of  C2*,  570;  C3+,  713;  C/,  605;  642;  ,  496;  and  --111  kJ 

mol'^,  consistent  with  higher  stability  for  the  odd-n-value 
cation  clusters. 

NR  Energy  Deposition.  The  differences  between  these  AHf(Cn^) 
values  and  the  previous  AHf(Cn°)  values  yield  IE(C„°)  values 
(underlined,  Table  I)  that  can  be  checked  further  against  the 
extent  of  precursor  dissociation  in  the  ’  X  spectra.  The 
appropriate  minimum  internal  energy  of  C„°  after 
neutralization  is  dependent  on  IE(Cn°)  -  IE  (target)  ,  with  target 
IE(Mg,Na,K)  =  6.1,  5.1,  and  4.3  eV.  As  an  example,®® 
neutralization  of  CH/  (IE  =  12.5  eV)  with  K  should  produce  CH*® 
with  E*oin  =  12.5  -  4.3  =  8.2  eV.  The  product  yields  are  CH^*,  0%; 
CH3'^,  53%;  and  CH2'^,  and  C^;  47%.  This  results  from 

dissociation  thresholds  of  4.6  eV  for  CH/  and  9.3  eV  for  CH3^,®^ 
1.1  eV  above  E*„in . 

The  neutralization  of  C3*  by  Na  (K)  should  produce  €3°  with 
E*niin  =  6.3  (7.1)  eV;  this  causes  53%  (75%)  Cj°  dissociation  that 
requires  7.5  eV,  or  1.2  (0.4)  eV  above  E*B,in.  Neutralization  of 
by  Mg  (Na)  should  produce  E\^„(C^°)  =  4.6  (5.6)  eV;  this  causes 
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59%  (94%)  C4°  dissociation  that  requires  5.4  eV,  or  0.8  eV  above 
(0.2  eV  below)  Neutralization  of  C5*  by  Mg  (Na)  should 

produce  E*ni„(C5°)  =  4.6  (5.6)  eV;  this  causes  28%  (69%)  €5° 
dissociation  that  requires  5.8  eV,  or  1.4  (0.4)  eV  above  E*B,in* 
Neutralization  of  Cg*  by  Mg  (Na)  should  produce  E*„in(C6°)  =3.7 
(4.7)  eV;  this  causes  75%  (96%)  dissociation  that  requires  3.6  eV 
or  0.1  (1.1)  eV  below  E*nin.  The  agreement  of  these  data  with 
each  other,  and  with  that  of  CH*,  is  -±0.3eV.  An  exce  1  is 
Cj®,  which  is  only  32%  dissociated.  The  excited  electronic 
states  of  the  small  Cj®  molecule  should  have  well  separated 
intersystem  crossings,  and  thus  some  could  have  sufficiently  long 
lifetimes  (>0.5  /ns)  to  reach  the  reionization  cell  before 
dissociation. 

Values.  Such  an  evaluation  for  the  negative  ions 
can  be  made  using  the  charge  exchange  ^NR"  spectra  of  Figure  3 . 
With  the  underlined  values  of  Table  I,  all  dissociations  are 
predicted  to  be  highly  endothermic  (850  kJ  mol'^  for  C2*)  , 
consistent  with  Figure  3.  From  the  formation  of  Cj'  is 
predicted  to  be  favored  by  198  kJ  mol'^.  From  C*',  differential 
enthalpy  values  are  ,  58;  C{ ,  0;  and  Ci',  124  kJ  mol'^.  From  Ci ^ 
,  these  values  are  C^',  96;  C3';  132;  C{ ,  0;  and  Ci’,  331  kJ  mol‘^. 
From  Cg’,  these  values  are  C5',  68;  C*',  38;  ,  0;  C^' ,  75;  and 

,  216  kJ  mol"^.  Eyler  has  recently  reported^®  electron  affinities 
for  C4®,  Cj®,  and  Cg®  of  1.3  -  1.5  eV,  far  lower  than  those  derived 
from  the  Table  I  values.  Note  that  dissociation  channels  other 


than  those  accompanied  with  C3  neutral  loss  are  also  prominent, 
contrary  to  the  photofragmentation  results.^®  These  AHf(C„  ) 
values  predict  minimum  energies  of  dissociation  of  ,  850;  C3', 
594;  C*',  630;  Cj",  526;  and  Cg*,  551  kJ  mol'S  consistent  with 
higher  stability  for  the  even-n-value  anion  clusters. 

Conclusions 

Although  multiple  isomeric  €3°  and  species  are  predicted 
to  be  stable,  those  prepared  from  a  wide  variety  of  linear, 
branched,  and  cyclic  (and  isotopically  labeled)  precursors  give 
NR  spectra  indicative  of  only  one  isomer,  presumably  the  linear 
one.  Consistent  with  results  from  gas-phase  ion 
chromatography,^®  CAD  cross  sections  indicate  that  C5*  and  Cg^  are 
also  linear,  but  that  C;"^  is  cyclic.  Values  of  AHf(Cn°,  Cn*,  Cn') 
and  IE(C„°)  selected  from  previous  experimental  and  theoretical 
studies,  plus  from  calculations  here,  can  be  fit  to  reaction 
enthalpies  consistent  with  these  NR  and  CAD  spectra,  indicating 
values  that  should  be  of  improved  relative  accuracy. 
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Table  I.  Heats  of  Formation  (298  K,  kJ  mol'^)  and  Ionization  Energies 


(eV) 

;  chosen  values 

underlined. 

Soecies 

iHTTcJi 

IE 

ThTTcTT” 

Cl  expt*’ 

717 

11.3 

I8O3 

595 

theory** 

717*= 

11.0 

1778 

theory* 

717c 

11.1 

1788 

C2  exph**'* 

832 

12.1 

1998 

505 

theory** 

815* 

12.1 

1983 

497 

theory** 

782 

12.1 

1950® 

462 

C3  expt*" 

837 

12.1 

2004 

649 

expt** 

13.0 

theory 

801^ 

11.4** 

1872* 

613 

theory® 

773 

10.4 

1776 

585 

C4  expt** 

971 

12. 6*’'*’ 

2187 

614 

expt* 

10.7 

2004 

theory** 

1022 

10.5 

2033 

665 

theory* 

1012 

11.3 

2102 

655 

C5  expt 

979* 

11.5* 

2024 

709 

expt 

12.3*= 

2012* 

theory** 

1030 

J.0.7 

2066 

760 

Cg  expt 

>1180** 

9.6* 

>2106* 

>784 

expt** 

9.7 

theory** 

1203 

9.8 

2149 

807 

C-C7  expt 

>1146“ 

>1923*^ 

>847 

theory** 

<1263" 

<2073" 

<964 

n-C7  expt 

>>1146*''" 

»2111"'*’ 

»847 

n-C7  expt 
theory** 


1263" 
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2227P 


964 


Footnotes  for  Table  I. 

-Using  aH£(C„°)  plus  vertical  electron  affinity  (EA)  values  from  ref. 

19  except  EA(Ci)  from  ref.  27;  EA(C3-C6)  from  ref.  12  are  lower  by  28- 
40  kJ  mol'^.  -From  ref.  27.  -Reference  value  used  in  the  theoretical 
calculations.  -Calculated  from  data  (total  electronic  energy, 
fundamental  frequencies,  IE,  etc.)  of  ref.  lla,  with  thermal 
corrections  and  heat  of  atomization  (binding  energy,  0  K)  scaled  by 

I. 1.  -As  in  d  from  our  calculations  at  the  MP4  (STDQ)  level;  IE(C2°) 
from  the  total  energy  of  the  C2*  quartet  state,  corrected  by  the 
excitation  energy  E('‘ng)  -  E('‘2g‘)  =  1.2  eV  from  the  MCSF-SCEP 
calculations,  ref.  42;  calculated  for  the  rhombus  structure. 
-Experimental  value  of  815  kJ  mol'^  from  Urdahl,  R.S.;  Bao,  Y.; 

Jackson,  W.M.  Chem .  Phvs .  Lett .  1991,  178 .  425-428.  *From  the  present 
experimental  data  and  the  underlined  IE  values.  l|From  ref.  28.  iFrom 
ref.  18c.  -^Chase,  M.  W. ;  Davies,  C.  A.;  Downery,  J.  R. ;  Frurip,  D. 

J. ;  McDonald,  R.  A.;  Syverud,  A.  N.  J.  Phvs.  Chem.  Ref.  Data  1985,  14, 
Supplement  1:  JANAF  Thermochemical  Tables,  Third  Edition.  -Drowart, 
J.;  Burns,  R.  P. ;  Demaria,  G. ;  Inghram,  M.  G.  J .  Chem .  Phvs .  1959,  31 . 
1131-1132.  -Linear  isomer  (ref.  lla).  -Assuming  the  cyclic  isomer  to 
be  the  more  stable.  -Using  the  theoretical  IE  value  from  ref.  18c. 
SRef.  28;  lE(linear)  -  lE(cyclic)  =  1.6  eV  from  theory  and  1.9  eV  from 
experiment,  ^sing  the  IE  value  f rom  Koopman ' s  theorem. 


Table  II  .  Calculated  total  energies  of  carbon  cluster  ions  and  neutrals. 


total  energy* 


Ion/ 

Neutral 

<S^>^ 

SCF  MP2 

MP4(SDQ)  MP4(SDTQ) 

ZPVE^ 

C 

('p) 

0 

-37.598730  -37.669305 

-37.694076  -37.694885 

C 

(^P) 

2.01 

-37.689049  -37.745023 

-37.763780  -37.764302 

C2 

( V) 

0 

-75.393686  -75.718201 

-75.720480  -75.757361 

11.4 

C2 

2.01 

-75.465389  -75.665165 

-75.681075  -75.691228 

14.3 

C3 

('V) 

0 

-113.378370  -113.726307 

-113.743722  -113.776683 

23.2 

C3 

‘X) 

2.21 

-113.283383  -113.645172 

-113.659625  -113.697914 

20.0 

c« 

('V) 

0 

-151.166223  -151.758917 

-151.628814  -151.859454 

41.9 

C, 

(X) 

2.28 

-151.214780  -151.618959 

-151.620997  -151.677634 

42.4 

C* 

(X) 

0 

-151.175999  -151.659151 

-151.673560  -151.708357 

35.3 

c*- 

(2p) 

0.76 

-37.291796  -37.336540 

-37.357692  -37.358019 

c*- 

(*P) 

3.75 

-37.158039  -37.167617 

-37.170724  -37.170734 

0.80 

-74.922290  -75.161871 

-75.170107  -75.188772 

13.  S 

Cz" 

(X) 

3.75 

-75.093854  -75.242087 

-75.262514  -75.268205 

13.2 

C3" 

(X)* 

0.85 

-112.938031  -113.301503 

-113.311787  -113.343776 

23.4 

C3" 

(X)*" 

0.85 

-112.936514  -113.303785 

-113.310393  -113.350904 

C3* 

(X")' 

0.78 

-112.910338  -113.321953 

-113.315460  -113.383621 

22.5 

C3" 

0.78 

-112.908702  -113.323758 

-113.317587  -113.394710 

C3" 

3.99 

-112.977334  -113.231845 

-113.255894  -113.272173 

23.1 

C/ 

(X*) 

0.89 

-150.828454  -151.237934 

-151.280966  -151.307179 

32.2 

C," 

(%) 

4.14 

-150.776199  -151.237982 

-151.264185  -151.308558 

36.2 

C/ 

(^A2) 

2.67 

-150.814654  -151.181499 

-151.219405  -151.240990 

35.1 

C/- 

(*A") 

3.86 

-150.831167  -151.246165 

-151.265529  -151.293238 

36.1 

1.37 


Footnoces  for  Table  II  . 

•Geometry  optimizations  and  energy  calculations  with  the  6-311G(d)  basis  set,  energy  in 
hartrees . 

^Spin  contamination  value  (ref.  40). 

'Zero  point  vibrational  energy  from  6-311G(d)  harmonic  frequencies,  unsealed,  kJ  raol"^. 
^Geometries  optimized  at  the  UMP2/6-311G(4)  level. 

•First-order  saddle  point,  single  imaginary  frequency. 

^Second- order  saddle  point,  two  degenerate  imaginary  frequencies. 

’Quartet,  state  unassigned. 


Table  III.  ■^NR'^  Spectra  of  , 


Ca(~85%T)  V02(70%T)  . 


Precursor 

m/z  12 

24 

36 

48 

n-C^He 

5 

16 

56 

23 

n-c^cie 

6 

16 

60 

19 

C" C^CIq 

6 

14 

60 

21 

C“CgHg 

5 

19 

56 

19 

C— CgClg 

7 

17 

58 

18 

“Transmittance  from  eso’C  Ca  vaporization. 
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Table  IV.  Cross  Sections  for  CAD  and  NR  Spectra. 


Cluster 

p  +  a 

^  +  b 

+  C 

^  0  d 

C„'  “ 

-♦c  ^ 

^<n 

■*C<n' 

Cl 

- 

15.4 

5.2 

- 

C2 

0.28 

12.6 

0.90 

2.1 

<1 

C3 

0.42 

9.8 

0.23 

1.8 

18 

C, 

1.8 

7.7 

0.24 

0.90 

25 

C5 

1.8 

6.3 

0.09 

0.60 

37 

Cg 

2.8 

4.9 

0.14 

0.30 

20 

C7 

0.80 

— 

— 

— 

— 

*C„*  CAD,  02(50%T)  ,  efficiency:  total  C^"^  product  ion  abundance,  in  % 
of  unattenuated  precursor  ion  abundance  (±20%  relative) .  -C„* 

neutralization,  Xe(70%T)  ,  efficiency:  total  C<n°  flux,  %  of 
unattenuated  precursor  ion  abundance  (±10%  relative)  .  -C^'^  charge 

reversal,  C6H6(30%T),  efficiency:  total  C^^'  abundances,  %  of 
unattenuated  precursor  abundance  (±20%  relative).  -C„°  [from  Xe(70%T) 
neutralization  of  C^^]  reionization,  02(70%T),  efficiency;  total 
abundances,  %  of  total  C<j,°  flux  (±20%  relative)  .  dissociation 

efficiency  in  charge  reversal  spectra,  C6H6(30%T)  ;  total  C<n'  product 
ion  abundance,  %  of  survival  Cn". 


Figure  Legends 


Figure  l 


Figure  2 


Figure  3 


Figure  4 


Calculated  geometries  of  C*'^:  i!/  »*  ^A"(D2h  -  C.)  ; 

2l,  *B,(D«h  -  Cah)  ;  4*,  ^AaCCzv);  Sl,  Bond  lengths  are 

in  Angstroms. 

CAD  spectra,  02(50%T)  ,  identical  (±10%)  for  the  same  ion 
from  the  listed  precursors:  (A)  €2^  from  n-CgHg,  n-C^Clg,*  (B) 
C2*  from  n-CaHg,  n-C^Clg,  c-CsHg,  c-CsClg;  (C)  C^*  from  n-C^Hg, 
n-C^H^Bra,  n-C^Clg,  i--C4H8,  i-C^HaBr,  i-C^HaBrj,  C-CgClg,  c-CgHg, 
and  c-CgClg;  (D)  C5*  from  c-CgClg;  (E)  Cg^  from  c-CgClg  and  c- 
CgBrg,*  and  (F)  C/  from  toluene  and  cycloheptatriene. 

Charge  reversal  *NR'  spectra,  C6H6(30%T),  identical  (±20%) 
for  the  same  ion  from  the  listed  precursors;  (A-F)  Cj^-Cg"^, 
respectively,  as  in  Figure  2,  except  C/  is  from  n-C^Clg,  i- 
C^HaBr,  and  s-CgClg.  Higher  trace  is  amplified  tenfold. 

spectra,  Y/02(70%T),  identical  (±10%)  for  the  same  ion 
from  the  listed  precursors;  (A)  Ci^,  (B,G,L)  C2*,  (C,H,M) 

C3%  (D,J,N)  C4\  (E,J,0)  C5%  (F,K,P)  with  (A-F)  Y  = 

Xe(70%T),  (G-K)  Y  =  Na(85%T),  (L-P)  Y  =  K(90%T) .  Precursors 

as  in  Figure  2,  except  C3*  is  from  n-CjHg  and  n-C^Clg  and  also 
from  c-CgHg  and  c-CjClg  with  K  neutralization;  C*^  is  from  n- 
C^Clg,  i-C4H7Br,  and  c-CgClg  and  also  from  c-CgClg  with  Na  and 
K  neutralization.  Numeric  values  for  precursors  are  the  % 
of  the  value  with  no  collisions  (value  for  Mg  neutralization 
in  parenthesis.)  Product  values  are  relative  peak  areas, 
corrected  for  dissociation  after  reionization  indicated  in 
the  Xe  spectra,  averaged  for  the  multiple  measurements. 


